Interstitials, Vacancies and Dislocations in Flux-Line Lattices: A Theory of Vortex 

Crystals, Supersolids and Liquids 
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We study a three dimensional Abrikosov vortex lattice in the presence of an equilibrium con- 
centration of vacancy, interstitial and dislocation loops. Vacancies and interstitials renormalize the 
long-wavelength bulk and tilt elastic moduli. Dislocation loops lead to the vanishing of the long- 
wavelength shear modulus. The coupling to vacancies and interstitials - which are always present 
in the liquid state - allows dislocations to relax stresses by climbing out of their glide plane. Sur- 
prisingly, this mechanism does not yield any further independent renormalization of the tilt and 
compressional moduli at long wavelengths. The long wavelength properties of the resulting state 
are formally identical to that of the "flux-line hexatic" that is a candidate "normal" hexatically 
ordered vortex liquid state. 



I. INTRODUCTION 

Both disorder and thermal fluctuations strongly affect 
the properties of the vortex array induced in terpe II su- 
perconductors by an external magnetic field.atx One 
of the most striking consequences of thermal fluctua- 
tions, particularly pronounced in the high-T c materials, 
is the resistive vortex liquid state,u located between the 
H C 2 (T) line and the vortex solid in_the magnetic field (H)- 
temperature (T) phase diagram. Q Upon field cooling, a 
vortex liquid freezes into an Abrikosov vortex solid. The 
nature of the freezing transition and of the resulting vor- 
tex solid phase depends on the amount of disorder present 
in the material. In dirty samples the vortex solid has been 
described as a "vortex glass" J3 and its translational cor- 
relation length is limited by disorder to a finite value.Q In 
three dimensions, the low temperature vortex glass solid 
is expected to be a true superconductor with a vanishing 
linear resistivity. For weak disorder, the vortex solid state 
is expected to be a topologically ordered "Bragg" glass 
state in three dimensions, with logarithmically gnawing 
vortex displacements, but bound dislocation loopsJj The 
freezing transition-of the vortex array has been observed 
to be first orderoO in ultra .clean samples and continuous 
in dirty superconductors.B~ll3 

In very clean samples, where the disorder-limited 
translational correlation length is thousands of intervor- 
tex lattice constants, the low temperature phase can be 
well approximated by a vortex lattice. Within an elas- 
tic description, the primary low temperature excitations 
of the vortex lattice are phonons, characterizing two- 
dimensional displacements of vortex lines from their pre- 
ferred lattice positions. As the temperature (or field) is 
raised towards the melting transition, other excitations 
become important. By definition, these are defects in 
the vortex lattice, i.e. they are not describable in terms 



of single-valued vortex displacements. These line defects 
are dislocations, disclinations, vacancies and interstitials 
and must be included in the model for a complete de- 
scription of the melting of the vortex solid and of the 
properties of the resulting vortex liquid state. 

At low temperatures, in a well-ordered Abrikosov lat- 
tice statejj these defects are bound, as the energy of an 
isolated line defect diverges with system size. At higher 
temperatures entropy can, however, drive a prolifera- 
tion of these line defects,, i n . analogy with the melting 
of two-dimensional solids .olij Two melting scenarios are 
possibleJla (i) a one-stage first order transition from a 
solid to an isotropic vortex liquid where both dislocations 
and disclinations unbind simultaneously (as it occurs in 
the melting of ordinary 3d-splids) , (ii) a two-stage, pos- 
sibly continuous transitionlij, where dislocations unbind 
first, leading to a hexatic flux- line liquid with residual 
bond-orientational order, vanishing shear modulus, but 
finite hexatic stiffness. This first transition would then be 
followed by a proliferation of disclination loops, thereby 
completing the transition into an isotropic vortex liquid. 
The first stage of this second scenario for the melting of 
the Abrikosay lattice was first suggested by Marchetti 
and Nelson.EJ While avoiding the subtle question of the 
melting transition itself, they adapted_the method de- 
veloped long ago by Nelson and TonerllZl to describe the 
vortex line-liquid state. Marchetti and Nelson described 
the hexatic vortex liquid as a vortex lattice with an equi- 
librium concentration of dislocation loops, treating the 
latter in the Debye-Huckle approximation. Through de- 
tailed calculations, they demonstrated that dislocations 
drive the long wavelength shear modulus of the system 
to zero and computed the effective hexatic stiffness of the 
resulting orientationally ordered vortex liquid.t^l 

Vacancies and interstitials constitute another class of 
defects that play an important role in solids. In ordinary 
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crystals and in 2d vortex lattices, these point defects cost 
finite energy and are therefore present in finite density, 
at any finite temperature. While their static effects in 
these systems are minimal, their density represents an 
important hydrodynamic mode which must be included 
in the correct description of crystal hydrodynamics S3 

In strong contrast, in vortex lattices, vacancies and 
interstitials are line defects with energy proportional 
to their length and thereby diverging in bulk samples. 
Hence one expects that at low temperatures their average 
density vanishes. At higher temperatures, this positive 
energetic contribution to the free energy can, however, 
be compensated by a negative entropic contribution as- 
sociated with line wandering, which also scales with the 
defect length, in analogy with the Kosterlitz-Thouless 
pictureJi-3 These considerations allow two thermodynam- 
ically distinct crystalline phases, with a sharp phase tran- 
sition between them. While in both phases dislocations 
and disclinations are bound, and as a consequence there 
is long-range translational order (true Bragg spots in an 
X-ray scattering experiment) and a finite shear mod- 
ulus, line vacancies and interstitials are bound in the 
low temperature crystal, but have proliferated in the 
high temperature crystal phase. A thermodynamically 
sharp distinction between these two crystal phases in 
three-dimensional vortex systems was first emphasized by 
M.P.A. Fisher and LeetLa based on the mathematical cor- 
respondence between vortex lines and world-lines of two- 
dimensional quantum bosons. In this mapping the low 
temperature vortex crystal maps onto a 2d Wigner crys- 
tal and the high temperature vortex solid corresponds to 
the quantum supersolid phase of 2d bosons, with vacan- 
cies and interstitials in its ground state. The quantum 
supersolid is quite exotic, in that it is both crystalline 
and superfluid. Correspondingly, due to the finite den- 
sity of vacancies and interstitials in the vortex supersolid, 
vortex lines can move arbitrarily far and entangle, as in 
a vortex liquid, and therefore this phase exhibits finite 
linear resistivity. 

While experiments seem to rule out the existence of 
an equilibrated vortex supersolid phase in bulk 2d quan- 
tum crystals, based on detailed calculations, Frey, Nelson 
and D.S. Fisher ,E3 have argued that such a phase is more 
likely to exist in flux-line arrays at high fields because of 
the layered structure of high temperature superconduc- 
tors. These authors conclude that a vortex supersolid 
phase will certainly exist in anisotropic superconductors 
for magnetic fields above a decoupling field B x where 
vortices in different Cu02 layers are essentially decou- 
pled by thermal fluctuations E3 Furthermore, even if an 
equilibrium supersolid phase was absent, an appreciable 
noncquilibrium density of vacancies and interstitials may 
still be present in a flux-line lattice, when the vortex ar- 
ray undergoes a first order freezing transition upon cool- 
ing in a constant field. 

In this paper we study the effects of vacancies and in- 
terstitials within the vortex supersolid and liquid phases. 
As discussed above, in the supersolid phase, aside from 



being responsible for its existence, these defects provide 
a mechanism for vortex line wandering and consequently 
for its finite resistivity. They also are important degrees 
of freedom, in addition to phonons, that must be incor- 
porated in the correct description of the supersolid. In 
Sec. Ill we construct a model of a vortex supersolid, as an 
clastic lattice with an equilibrium concentration of un- 
bound vacancies and interstitials. We compute the flux- 
line density correlation functions, that characterize the 
static equilibrium properties of this phase, and extract 
from these the effective elastic moduli of the supersolid 
phase. We find that the long-wavelength shear modulus 
is unaffected by the fluctuations in the density of vacan- 
cies and interstitials. This result is consistent with the 
vacancy-interstitials' inability to relax shear, and con- 
firms the finiteness of the supersolid shear rigidity which 
distinguishes it from a vortex liquid. We also compute 
a finite, downward renormalization of the compressional 
and tilt moduli by vacancy and interstitial density fluc- 
tuations. We demonstrate that as a consequence of the 
reduction of the effective tilt modulus, the flux-line wan- 
dering is enhanced and, analogously to a vortex liquid, 
a bulk vortex supersolid is always entangled, i.e., it does 
not exhibit longitudinal superconductivity. 

The existence of a vortex supersolid allows for two sce- 
narios for the melting transition into the vortex liquid 
state. At low magnetic fields, we expected this transition 
to be directly from the low temperature, non-supersolid 
crystal into a vortex liquid phase. As discussed above, 
the vortex liquid state itself can be either a fully disor- 
dered isotropic liquid or a bond-orientationally ordered 
liquid that can further disorder into an isotropic liquid 
via disclination unbinding. 




Fig. 1. Schematic B — T phase diagram illustrating the flux 
line crystal, vortex "supersolid" (guaranteed to exist at 
fields much higher than the decoupling field B x ), the 
orientationally ordered hexatic vortex liquid, and the 
isotropic vortex liquid. 

Alternatively, the non-supersolid crystal can first un- 
dergo a transition into a vortex supersolid by a prolif- 
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eration of vacancies andjaterstitialsj^l and subsequently 
melt into a vortex liquidE3 These two scenarios have-been 
dubbed as type I and type II melting, respectively,c2l and 
are illustrated in Fig.l. 

Very recent experiments investigating vortex penetra- 
tion through surface barriers in the presence of a trans- 
port current in clean BSCCO samples have indeed indi- 
cated that the nature of the melting line may change at 
high field, even in the absence of point disorder In fact 
these experiments suggest the existence of an intermedi- 
ate phase between the solid and the liquid at high fields. 
This observed regime lies above the first order melting 
line determined by equilibrium magnetization measure- 
ments, but apparently exhibits a certain amount oLsjx- 
fold periodicity, as indicated by neutron scattering. Ell It 
is therefore a candidate for either the supersolid or the 
hexatic liquid phases studied here. 

Although the nature of the type I .and type II melting 
transitions should be quite different ,E3 in either case va- 
cancies and interstitials should proliferate in the resulting 
vortex liquid state.tJ Hence, in addition to dislocations, 
vacancy and interstitial defects must be taken into ac- 
count for the proper description of a vortex liquid. They 
were not, however, explicitly included in the calculations 
of Ref. [L6|. When vacancies and interstitials are absent, 
dislocation loops in the vortex line-lattice are restricted 
to lie in the plane defined by the magnetic field axis and 
their Burger's vector, and can only glide (see Fig. 5), but 
not climb. One of the consequences of this is that the 
effective tilt modulus of the vortex liquid is not renor- 
malized by dislocations and in the model of Ref. [l6| is 
therefore identical to that of the vortex lattice. 



The results presented in Sec. IV remedy the limitations 
of Ref. |l^, by explicitly including vacancies and intersti- 
tials in the description of a vortex liquid. Vacancies and 



interstitials renormalize the tilt and compressional mod- 
uli and allow dislocation loops to climb out of their glide 
plane by absorbing and emitting these defects. At long 
wavelengths, however, the coupling of dislocations to va- 
cancies and interstitials does not yield any additional in- 
dependent renormalization of the tilt and compressional 
moduli. Whether this is a general property of the vor- 
tex lattice, or an artifact of the quadratic model and 
the Debye-Huckel approximation used here, remains an 
open question. Our work yields a complete description 
of a bond-orientationally ordered hexatic vortex liquid, 
whose tilt and bulk moduli are renormalized by disloca- 
tions, vacancies and interstitials, and whose shear modu- 
lus is driven to zero by the proliferation of these defects. 

Our results characterizing the properties of the vari- 
ous phases shown in Fig.l are summarized in the two 
tables below. We stress that hexatic I and hexatic II are 
not two distinct hexatic phases, but rather two differ- 
ent regimes within the same hexatic phase, distinguished 
by a high and low concentration of vacancy/interstitial 
defects, respectively. The crossover between these two 
regimes is indicated in Fig.l by a dotted line. The re- 
sults quoted in Table 1 for the hexatic I phase are those 
obtained by Marchetti and Nelson in Ref. |l6| assuming 
that no vacancy/interstitial defects are present in the 
system. These results, however, only apply well below 
the dotted line in Fig.l, very close to the solid-hexatic 
phase boundary. As the transition to the isotropic liquid 
is approached, the large number of interacting disloca- 
tion loops present in the hexatic will invitably generate 
vacancy /interstitial defects as well, leading to the break- 
down of the model of Ref. 16. The type of long-range 



order present or absent in each of these phases is sum- 
marized qualitatively in Table 2. 
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Table 1. A summary of our results for the elastic constants and the winding number in the various phases of Fig. 1. 
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Translational 


Orientational 


Longitudinal 




LRO 


LRO 


Superconductivity 


crystal 


yes 


yes 


yes 


supersolid 


yes 


yes 


no* 


hex I 


no 


yes 


no 


hex II 


no 


yes 


no 


liquid 


no 


no 


no 



Table 2. This table displays the presence or absence of translational and orientational long-range order (LRO) as 
well as longitudinal superconductivity in each phase. The asterisk serves to emphasize that although the supersolid 
does not exhibit longitudinal superconductivity, the degree of screening of longitudinal currents in this phase is 
substantially different from that of the vortex liquid phase (see Table 1). 



II. ELASTIC PROPERTIES OF DEFECT-FREE 
VORTEX LATTICES 

A. Model 

We begin by recalling the properties of an ordered, 
defect-free vortex lattice, which we expect to be stable at 
low temperature and field in a clean superconductor .0 As 
discussed in the Introduction, long-scale degrees of free- 
dom of this system are uniquely characterized by a single- 
valued vortex displacement field, u(r). With the convec- 
tion in which the external magnetic field Ho (aligned 
with the c axis of the superconductor) points along the 
z axis, the two-dimensional vector displacement u(rj_, z) 
is confined to the xy plane. 

The long-wavelength properties of a triangular flux- 
line lattice are characterized by the elastic free energy 
functional 



i*latt = ^ / </l' 



where 



2c 66 w 4 2 + (en - 2c 66 )u 2 kk + c 44 ((9 z u) 



(2.1) 



(2.2) 



is the symmetrized two-dimensional strain tensor, with 



d a Uj 



(2.3) 



a 3 x 2 hybrid strain tensor. Greek indices take on the full 
three-dimensional set of labels x, y, z, and Latin indices 
are reserved for the purely two-dimensional set x,y. We 
will use this notation throughout the manuscript. The 



parameters cqq, c\\ and C44 appearing in the Eq.2.1 



the shear, compressional and tilt modulus, respectively. 
In contrast to ordinary crystals, in a flux-line lattice 
vortex interactions extend over a range of order of the 
London penetration depth, A, which can be quite large, 
especially in high T c superconductors. As_emphasized 
in the extensive literature on the subject,Era on scales 
shorter than A, this leads to wavevector-dependent elas- 
tic moduli. For a detailed comparison with experiments, 



inclusion of these nonlocal elastic effects can be impor- 
tant, especially at high fields, and they can be eas- 
ily incorporated into our results by replacing all of the 
bare elastic moduli by the proper wavevector-dependent 
expressions. E£l 

When the lattice contains no vacancies nor interstitials, 
the number of flux lines equals the number N of sites in 
the triangular lattice. On the average, the flux lines are 
aligned with the external field and the equilibrium mag- 
netic flux density field is given by Bo = zBq = zc/jqUq, 
where <fio = hc/2e is the flux quantum, no = N/A = l/a c 
is the equilibrium number density of vortex lines and a c 
the area of the primitive unit cell. 

Fluctuations in the local induction <5B(r) = B(r) — Bo 
can be described in terms of fluctuations in the areal den- 
sity of flux lines, <5n(rj_, z) — n(r±, z) — no, and of a tilt 
vector field t(rx, z), with the relation 

6B Z = o <5n, (2.4) 
B± = o t, (2.5) 

valid in the long wavelength qX << 1 limit .0 In the ab- 
sence of vacancies and interstitials, the areal density of 
flux lines and their orientation relative to the applied field 
direction are entirely determined by the local strains, ac- 
cording to 



Sn/n = -8A/A = 
t/n = d z u. 



(2.6) 
(2.7) 



The condition V • B = translates into a "continuity" 
constraint for the flux lines, 



d z Sr 



V_l - t = 0. 



(2.8) 



As can be seen from Eqs.2.6 and 2.7, this continuity con- 



straint is identically satisfied in the defect-free vortex lat- 
tice, where the displacement u is single- valued. 



B. Correlation and response functions 

Thermal fluctuations in the density and tilt field are 
characterized by the density-density correlation function 
(the structure factor), 
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S{q±,q z ) = i(<5n(q)5n(-q)), 

and the tilt field correlation function, 
1 



Tij(q±,q z 



where 



V 



(ti(q)ti(-q)) 



(2.9) 



(2.10) 



T L (q± , q z )P l L j (q± ) + T T (q ± , q z )P7 (q ± ) , 



(q±) 



f«(qx) = % 



(2.11) 
(2.12) 



are the longitudinal and the transverse projection oper- 
ators, respectively, = qj_/g±, and V is the volume 
of the superconductor. In light of the constraint, Eq. 2 



the longitudinal part of the tangent field correlator is 
proportional to the structure function, with 



TL(q±,q z 



-§-S{q±,q z ). 

q± 



(2.13) 



The brackets (. . .) in above expressions and throughout 
the paper indicate a thermal average with a Boltzmann 
weight e~ F ' kBT /Z, with F the free energy functional and 
Z = r TYa,ce[e~ F / kBT ] the corresponding partition func- 
tion. In the defect-fr ee ( non-supersolid) vortex lattice, F 
is given by -Flatt, Eq.2.1, and we denote the correspond- 
ing thermal averages by (. ..)o- Using Eqs.2.£ and 2.7, 
the structure function and the tilt field correlation func- 
tion can be expressed in terms of thermal averages of the 
phonon field u, and are therefore easily computed, with 
the result, 



S°(q) = 



V 

n\k B Tq\ 



q± • u(q) 



cuqi 



and 



T°(q) = 



^J^(4 L )< Ui (q)u i (-q)) 



(2.14) 



(2.15) 



n^k B Tqi 



C669, 



c^qi 



The structure function S*(q) can be probed in a neu- 
tron scattering experiment. The tilt correlation function 
Ty(q) is directly connected to the experimentally mea- 
surable linear magnetic susceptibility tensor, Xiji ^^ ac ~ 
cording to 



^Xij(q±,qz) = + 



4irk B T 



Zk(?x s &). (2-16) 



Equation [2.16 holds in the long wavelength qX << 1 
limit. A more general relationship between susceptibility 
and tilt correlation function that applies at scales shorter 



than A in an anisotropic material can be found for in- 
stance inRcf. [55]. The first term on the right hand side 
of Eq.| .16| represents a perfect diamagnetic (negative) 
Meissner response, which, in a mixed state is consider- 
ably reduced by the "normal" paramagnetic vortex tilt 
response Ty, contained in the second term. The linear 
susceptibility relates the transverse flux density <5B^(q) 
induced by an external perturbation field <5Hj_(q), ap- 
plied perpendicular to the field H = zH responsible 
for the onset of the vortex state, with 



(2.17) 



From Eq.2.15, appropriate for a perfect, defect-free 



flux-line lattice, we observe that the long wavelength 
limit of the transverse part of the tilt field correlation 
function is nonanalytic, with 



T£(<Zi 



o,q z ) 

* = 0) 



nlk B T 



C44 







(2.18) 
(2.19) 



The nonanalyticity of the tilt correlation function re- 
flects a drastically different linear response of the defect- 
free vortex lattice to two types of transverse field 
perturbations £3 

The transverse field response in the limit q± < < q z — > 
0, corresponds to a tilt perturbation of the flux lines, 
induced by an applied transverse field Vdq ± _oSH±(q z ) 
that is spatially homogeneous in the xy plane. The 
corresponding long wavelength transverse susceptibility, 
XT = P^Xtj, is given by 



lim Xt (q± = 0,q z 



1 

47T 



B 2 



4irc. 



44 



(2.20) 



If the seco nd term in brackets on the right hand side 
of Eq.^20| were absent, the superconductor would ex- 
hibit perfect screening of the transverse perturbation. 
Such a behavior can for instance occur in flux-line ar- 
rays pinned by aligned damage tracks. In the pres- 
ence of such correlated disorder, the vortex lattice is 
replaced by a thermodynamically distinct "Bose" glass 
phaseJ23 that is characterized by a divergent tilt mod- 
ulus C44 and exhibits a transverse Meissner effect, with 
lim gi _>o Xt(?z) = — 1/47T. In the absence of anisotropic 
pinning of the vortex lattice, such a perfect transverse 
diamagnetic response is spoiled by the finite vortex tilt 
response, that leads to only a parti al screening of the 
transverse field, displayed in Eq. [2.20 

Conversely, the limit q z << q± — > 0, describes a mag- 
netic response to a transverse field VS qz _oSH±(c[±) that 
is homogeneous along the z-axis, but is spatially varying 
in the sy-plane. The induced z-directed screening cur- 
rents lead to a shear perturbation of the flux-line array, 
with the response in a defect-free lattice given by 



lim Xt(Q-1 



0) 



1 

47r' 



(2.21) 
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Thus the flux lattice exhibits perfect screening in re- 
sponse to this z-independent transverse perturbation. 
Since the screening currents involved in the shear pertur- 
bation run parallel to the applied field Ho || z, a perfect 
Meissner response to a shear perturbation has also been 
termed longitudinal superconductivity. It follows directly 
from the fact that, in contrast to a liquid, a vortex lattice 
is characterized by a finite shear modulus cqq. 

For comparison, we recall that in a flux-line liquid, the 
long wavelength limit of the transverse part of the tilt 
correlation function is analytic, with 



lim T l ' quid 

q±^0 1 



(Q±,Qz 



= 0) = lim Tip 

q^O 1 



liquid / 



0,1z 



nlk B T 

liquid 



(2.22) 



As expected, a vortex liquid, being qualitatively identical 
to the normal state, albeit highly conductive, exhibits 
neither transverse Meissner effect, nor longitudinal su- 
perconductivity. We will return to this point again in 
Sec.0. 

Finally, we note that in the latter q z << q± — ► limit, 
the transverse part of the tilt correlation function corre- 
sponds to the world- lines. winding number (W 2 ) studied 
by Pollock and CeperleyEa in their path integral approach 
to the superfluid transition in quantum boson systems. 
In such an approach, dating back to Feynman, the super- 
fluid phase is identified with an entangled state of boson 
world-line trajectories, 



lim T T (q x ,q z = 0) = (W 2 ) , 
q±^0 

hn s 
to ' 



(2.23) 
(2.24) 



and n s the boson superfluid density. The well-known 
identification of physical parameters under the boson- 
vortex mapping is given by 



h(3^ L , 



(2.25) 
(2.26) 
(2.27) 



with 777 the boson mass, t\ the core energy per unit of 
length of a single vortex line, (3 the inverse boson temper- 
ature and L the vortex sample thickness. Utilizing this 
mapping, together with the approximate local expression 
for the tilt modulus C44 = 77oei, we reassuringly find that 
the defect-free flux-line lattice (in which vortex lines do 
not entangle and the sample exhibits longitudinal super- 
conductivity) corresponds to the "normal" boson crystal 
with vanishing superfluid density n s — 0. The entangled 
flux-line liquid, on the other hand, corresponds to the su- 
perfluid phase of bosons with ti s = no, and does not ex- 
hibit longitudinal superconductivity. We stress, how ever, 



that the mapping described above, as well as Eqs. 2.23 



and 2.24, only apply to a model of the vortices that ne- 



the field (z) directionHa A more general approach to the 
derivation of vortex liquid hydrodynamics, directly from 
a "kinetic theory" of interacting flux lines was developed 
in Ref. M 



Hall 
sensors 



5H 



± 



5H 



5H, 



Fig. 2. Sketch of an idealized experimental setup to probe 
longitudinal superconductivity. 

Figure 2 shows a sketch of an idealized experimental 
setup that could be used to probe longitudinal supercon- 
ductivity. An external field H is applied along the axis 
[z in our coordinate system) of a cylindrical sample and a 
uniform vortex state with flux lines running along z is set 
up. A current / confined to a wire running along the axis 
of the cylinder and producing an additional azimuthal 
magnetic field <5Hj_(rj_) = 27z x rj_ /(r\c), which pro- 
vides the shear perturbation described above, can there- 
fore probe the q z = field response of the vortex state. 
The longitudinal superconductivity can then be studied 
by measuring the induced azimuthal component B^r^) 
of the local induction inside the sample. This may be 
possible by placing radially-directed row of Hall sensors 
at the top of the sample. In the defect-free crystal, which 
is a longitudinal superconductor, we expect B^ = ev- 
erywhere in the bulk of the sample, deeper than the pen- 
etration length A from the surface of the cylinder. In con- 
trast in the supersolid and hexatic phases there will be a 
nonvanishing azimuthal response B^ir^) everywhere in 
the sample. 



III. ELASTIC PROPERTIES OF SUPERSOLIDS 



gleets the nonlocality of the intervortex interaction along 
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A. Model 

As discussed in the Introduction, we expect that upon 
increasing the temperature and the external magnetic 
field, the defect-free vortex crystal will undergo a ther- 
modynamically sharp transition into a vortex supersolid, 
characterized by the coexistence of crystalline order and 
a finite—equilibrium density of vacancy and interstitial 
defectsj£9 Our goal here is to develop a continuum de- 
scription of the long wavelength clastic properties of such 
a supersolid phase. 

Once vacancy and interstitial line defects proliferate in 
the supersolid, their positions and orientations represent 
new and important low-energy degrees of freedom, inde- 
pendent of the lattice displacements characterized by the 
field u(r). At finite temperature, these defects will lead 
to fluctuations in the local induction, independent, but 
energetically coupled to, local elastic strains. To incorpo- 
rate such defect fluctuations, we adaptpthe hydrodynamic 
methods developed for vortex lines |3E!3 to the hydrody- 
namics of the vacancy and interstitial line-liquid. 

At finite density, i.e., within the vortex supersolid 
phase, the low-energy vacancy and interstitial configura- 
tions can be parameterized by z- directed conformations 
(r°(z),z), (with a = i,v denoting interstitial and va- 
cancy, respectively) as they traverse the sample along the 
direction of the applied field, in close analogy to vortex 
lines themselves. Long wavelength properties of a gas of 
N v vacancies and Ni interstitials can then be described 
in terms of a net areal density of defects, 



n d (r±,z) = rii(r±,z) - n v (r±,z) , 



and a two-dimensional tilt vector field, 



td(rj_,z) =ti(r_L,z)-t„(r_L,z) 



where 



i a (r ± ,z) = J2t (2) (r±-<(z)) 



(3.1) 



(3.2) 



(3.3) 



and 



N„ 



t a (r±,z) = 



dv a n {z) 
dz 



(3.4) 



Here nd{r±, z) represents the net number of defect lines 
crossing a unit area perpendicular to the field direction, 
while tdi(jj_, z) is the net number of defect lines crossing 
a unit area normal to the i-th direction, with i = x,y. 
Since the interaction among vacancy and interstitial de- 
fects is short-ranged, we expect that a liquid of such de- 
fects be characterized by a compressibility (inverse de- 
fect bulk modulus), Xi an d & finite tilt modulus, K. The 
corresponding long-wavelength free energy functional is 
therefore given by 



dr 



X -\6n d ) 2 + K{t d f 



(3.5) 



and nP, the mean net defect density 



with Srid = rid 
in equilibrium. 

In the presence of these defects, fluctuations in the lo- 
cal magnetic induction (or in the corresponding flux-line 
density n(r) and flux tangent field t(r)) can be brought 
about by changes in both the local lattice strains w a j 
and the defect densities rid and t d , as 



and 



Sn = —nowu + Srid 



ti = riQWzi + t dl 



(3.6) 



(3.7) 



In this case the continuity equation (2.8), arising from 
V • B = yields the nontrivial constraint 



d z 5n d + V_l • t d = 



(3.8) 



that defect lines cannot start nor sto p in side the sample. 
The elastic strain drops out from Eq.3.8, as it identically 
satisfies the constraint due to the single- valuedness of the 
displacement field u(r). 

Since motion of defects microscopically corresponds to 
hopping of vortex lines, we expect an energetic coupling 
between fluctuations in the density and orientation of 
defects and the elastic strain field. The lowest order cou- 
pling allowed by symmetry corresponds to the following 
interaction part of the free energy functional 



1 

no 



dr 



jSridV ■ u + j't d ■ d z u 



(3.9) 



where 7 and 7' are positive phenomcnological coupling 
constants with dimensions of elastic moduli. Ignoring a 
weak coupling to fluctuations in the local temperature, 
the free energy functiona ls F j nt and Fd, together with 
the elastic part -Fi a tt) Eq.2.1, completely determine the 
long-scale elastic properties of the vortex supersolid. 

It is important to note that the parameters 
en, C66, C44, 7, 7', x, and K, appearing in our roadel, are 
functions of the mean net defect density n^,L3 which, 
within the supersolid phase, can in principle be deter- 
mined through detailed microscopic calculations of the 
type presented in Ref. ^0|. At the vortex-crystal to su- 
persolid transition, we expect these parameters to dis- 
play a nonanalytic behavior as a function of the distance 
from the transition, |T — T xss \ (where T xss denotes the 
crystal-to-supersolid transition temperature), of the form 
illustrated for c 66 (n° d (T)) in Fig.3. 
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T 



xss 



T 



(3.11) 

In the above, the brackets denote a thermal average 
evaluated with the total supersolid free energy functional 
which includes both the elastic and the defects degrees 
of freedom, given by 



F SJ 



-Flatt + Fd + F i7 



(3.12) 



To ensure the condition of V • B = 0, these averages 
must be carried out under the nontrivial constraint that 
defect lin es do not start nor stop inside the sample, given 
in EqfJ. 

Utilizing this constraint to explicitly eliminate t^ in fa- 
vor of 5n d and re-expressing the strain tensor Wij in terms 
of the longitudinal and transverse single-valued lattice 
displacements 



Fig. 3. A sketch illustrating the behavior of the elastic 
shear modulus cqq near the vortex crystal to vortex su- 
persolid transition. The proliferation of vacancies and 
interstitials is responsible for the softening of the effective 
shear modulus inside the supersolid phase. 

In addition to these mean-field effects, the coupling 
of the elastic degrees of freedom to the fluctuations in 
defect density (around the average n d ) and orientation 
field, Eq.3.9, yields further renormalization of the elastic 



constants, study of which is in part the focus of our work 
here. 



B. Correlation and response functions 

The elastic properties of supersolids can be character- 
ized in a number of distinct ways, reflecting a variety 
of experimental probes that couple differently to the su- 
persolid degrees of freedom. The simplest of these, from 
the theoretical and experimental points of view, is the 
extension of the equal-time equilibrium cor relat ion f unc- 
tions S(q) and 7Y,-(q) 



defined by Eqs.gJ and |2.10| , re- 
spectively. They directly measure fluctuations in the lo- 
cal magnetic induction B(r), related to the fluctuations 
in the t otal vortex -line number and tilt densities, via re- 
lations (2.4) and (2.5). These latter quantities are de- 
termined by both the loca l ela stic stra in and the defect 
density fields through EqsJ3.' 



and 3.7, and lead to 



1 



s ss (q) = ^(lv„(q)-^(q)f) 



(3.10) 



Tff(q) = -((n w zl (q) + t dl (q))(n w zj (-q) + t dj (-q))) . 



UL(q) = q_L • u(q) 
it r (q) = (z x qjj ■ u(q) 



(3.13) 
(3.14) 



which are the independent finite wavevector elastic de- 
grees of freedom in the bulk, we obtain the total free 
energy characterizing the supersolid 



TP — 
± .s.s — 



d 3 q Jl 



(2tt) 3 



rT(q)| U T(q)| 2 + ^r i ( q )| Ui (q)| 



2n 2 
i 

"o 



79J 



" ' Qz - ) 5n d (q)u L (-q) 



K\td(q) 



QJ-J 



l'qztd{q)u T (-q) 



(3.15) 



In the above we have defined the transverse and longitu- 
dinal wavevector-dependent stiffnesses 



iMq) 

r £ (q) 



cnq± 



caq z 



(3.16) 
(3.17) 



and not surprisingly found that the transverse and the 
longitudinal degrees of freedom decouple. After re- 
expressing the S(q) and Ir(q) in Eqs.3.10 and 3.11 in 



terms of these same independent degrees of freedom, 
these correlation can be easily computed by inverting the 
corresponding 2x2 m atrice s that can be read off from the 
expression for F ss , Eq.3.15. For the structure function of 
the supersolid we thereby obtain 



S ss (q±,q z ) = y{\ -in qj_u L (q) +5n d (q)| 2 ) 



n k B Tq' 



(ciiqj + c AA g 2 z ) + (x^qj + Kq 2 z ) + 2( 7g j - jq 2 z ) 
{c liq \ + cmDix-^l + Kq 2 z ) - ( 7 ?1 ~ l'q 2 ) 2 



(3.18) 



Similarly, the transverse part of the tilt correlation function is given by 
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T} s {q±,q z ) = i(|tno?,UT(q) +t^(q)| 2 ) 



X +^in ^ c66+(? 2 (c44 _ 7 , w 



k B Tn 2 z . (3.19) 



It is convenient to define defect-renormalized elastic 
constants that characterize the effective elastic proper- 
ties of supersolids. One choice of such a definition, that 
naturally arises in experiments which measure fluctua- 
tions in the total magnetic induction B (e.g., neutron 
scattering experiments), is in te rms of th e static correla- 
tion functions given in Eqs. |3.18| and |3.19| , iden tifie d with 
their defect-free lattice counterparts, Eqs.2.14 and 2.15 



rijk B T 



S ss (q ± ,q z = 0) 



and 



nlk B T _ 



T^(q ± = 0,q z ) 



(3.20) 



(3.21) 



Similarly, the effective shear modulus c|| of the super- 
solid can be defined in terms of another equilibrium 
equal-time correlation function 



k B T 



■ 06 



(?x) 



= V ^(wijiqjWiji-q)) - (w u (q)w j3 (-q))] 
= V- l q\Pl (qx)<^(q) % (-q))| 9a =o , 



The typical behavior of the supersolid bulk modulus 
Cjj for en 7^ x _1 as a function of 7 is displayed in Fig. 4 
(C44 behaves similarly as function of 7'). The figure shows 
that the supersolid modulus grows (decreases) linearly 
with positive (negative) coupling, at weak coupling. At 
an intermediate positive value of 7, it reaches a maxi- 
mum at the value of the smallest of the moduli for the 
two subsystems. In the strong coupling regime (both 
positive and negative) the effective modulus decreases, 
vanishes and even changes sign, indicating an instability 
in the quadratic model of the supersolid. In our model, 
this instability is a signal of a (spurious) phase transi- 
tion within the supersolid phase, which in the case of the 
bulk modulus cf\ , corresponds to additional proliferation 
of defects and change in the lattice density. However, as 
we show in Appendix A, for a physical vortex supersolid, 
the elastic moduli (cg6, en, C44, x, K, 7, 7') that appear 
in our model are constrained to lie outside of this unsta- 
ble range. In the special degenerate case of en = \ , 

the effective bulk modulus of the supersolid behaves as 
r R _ 



qz=0 



(3.22) 



which can be easily evaluated by performing a Gaussian 
thermal average with the free energy f uncti onal, Eq. 3.15 



in the Boltzmann weight. Usi ng Eq s.BTj and 3.19 



gether with the definitions, Eqs. |3.20] and |3.21| , we find 



to- 



(cn + 7 )/2, with an analogous result for the tilt 



modulus C44. 

The calculation of c 6 
that, not surprisingly, vacancy and interstitial density 
fluctuations do not renormalize the shear modulus, and 



G 66' 



as defined in Eq.3.22, shows 



c«s(?-l) = c 6 6 



(3.25) 



In contrast to dislocations (to be considered in the next 
section), vacancies and interstitials are unable to relieve 
a shear stress. 



and 



cuX 



cn +27 



caaK - V 2 
c 44 + K - 2 7 ' 



(3.23) 



(3.24) 



From these expressions we observe that for vanishing cou- 
plings 7 = 7' = 0, the effective elastic moduli of the su- 
persolid are the corresponding moduli of the lattice and 
the liquid of defects, added in "parallel" — a physically 
appealing result. For weak coupling this implies that the 
effective bulk and tilt moduli of the supersolid are always 
smaller than or equal to the minimum of the correspond- 
ing moduli of the two subsystems, the lattice and the de- 
fects. Hence we find that for 7, 7' small compared to the 
elastic moduli of the lattice and the defects, fluctuations 
of vacancy and interstitial densities always reduces the 
effective longitudinal and tilt moduli of the supersolid, 
relative to that of the defect-free crystal. 




Fig. 4. The behavior of the supersolid bulk modulus as 
function of the coupling 7 between the elastic and the 
defect degrees of freedom, for c\\ = 1 < x^ 1 ■ 
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The presence of vacancy and interstitial defects also al- 
ters the response to a transverse magnetic field. The long 
wavelength limit of the transverse part of the tilt corre- 
lation function, T^ s (q±, q z ), remains nonanalytic, as in a 
defect-free crystal, and yields 



lim Xt (9-1 



lim n XT (q±>iz = o) = - - 



1 


1 - 


B 2 -I 


4-7T 


4?rcf 4 J 


1 


1 - 


B 2 i 


47T 





(3.26) 
(3.27) 



Since, as discussed above, these defects soften the lat- 
tice by decreasing both the longitudinal and tilt moduli 
the penetration of a transverse field is enhanced, as seen 
from Eqj3.26|. A more dramatic effect of the defects is 
the pres ence of the second term on the right hand side of 
Eq. |3.27 . In the vortex supersolid vacancy and interstitial 
defects allow flux-line wandering and entanglement, and 
as a result there is no perfect screening of weak trans- 
verse fields uniform along the z axis, and consequently 
the phase does not exhibit longitudinal superconductiv- 
ity. The supersolid is both crystalline (cqq ^ 0) and 
entangled {T^, s (q^,q z = 0) ^ 0), as argued by Frey, et 
al.E3 and consistent with nonvanishing winding correla- 
tion function 



(W 2 



lim T T {q±,q z = 0) 

q±^0 



(3.28) 
(3.29) 



We stress that the definitions of the effective elas- 
tic moduli characterizing a vortex supersolid is far from 
unique. This is related to the variety of experiments that 
probe different physical properties of the vortex super- 
solid. Instead of the corre latio n fun ctions S ss (q) and 
T^ s (q) defined in Eqs.3.10 and 3.11, we could have in- 
stead focused on 



5 ss (q) = ^(Ki(q) 



(3.30) 
(3.31) 



which probe only fluctuations in the lattice positions and 
not the total vortex density (related to magnetic induc- 
tion B). These latter correlation functions are more 
difficult to probe in a physical experiment, but can be 
straightforwardly measured in a numerical simulation. 
Making the identification between the effective elastic 
moduli of the supersolid and the co rresp ondi ng co rre- 



lation functions, in analogy with Eqs.3.2C and 3.21 
define 



n 2 k B T 



S ss (q ± ,q z = 0) 



and 



nlk B T _ - 



Ti s (q ± =0,q z ) 



(3.32) 



(3.33) 



Simple computation of >?22(q) and T^ s (q), together with 
these definitions leads toci] 



'-11 



c 44 
L 66 



Cll 
C44 
C66 



7 2 X > 
K ' 



(3.34) 
(3.35) 
(3.36) 



Another experimentally relevant way to probe elastic 
properties of vortex supersolids is through the linear re- 
sponse to a constant stress a a j applied at the boundaries 
of the system. For simplicity, we confine our discussion 
here to a two-dimensional stress cry applied to a bound- 
ary lying in the xy-pl&ne. To study the response, we first 
need to decide to which physical quantity does such stress 
couple. In a defect-free crystal, the answer is simple: the 
stress o~ij couples to the lattice strain My. In a supersolid, 
there is, however, a number of possibilities, depending 
on the nature of the experiment one seeks to describe 
(as was the case with the correlation functions discussed 
above). In a real (as opposed to a numerical) experiment 
the stress on the vortex lattice is produced through an 
electromagnetic interaction and therefore couples to the 
magnetic induction, which involves both the elastic lat- 
tice strain tensor and the defect contribution. Arguments 
similar to those found in Ref. |3^ indicate that the linear 
response of the supersolid can be studied by adding to 
the free energy F ss in Eq. |3.15 
stress, 



a part due to the external 



j dvcT,,in,, - -^—8 i: j8n d ) , (3.37) 



and then by minimizing it with respect to the indepen- 
dent elastic strain and defect degrees of freedom. 

Before proceeding with the calculation, some remarks 
are in order. Invoking the fluctuation-dissipation theo- 
rem, one might naively conclude that such a static re- 
sponse function is identical to the co rresp onding equal- 
time correlation function S ss (q) , Eq. 3. 10| that we stud- 
ied above. However, this is in fact not the case, in gen- 
eral, even for a defect-free crystal. To understand this 
difference, we consider for simplicity the case of a two- 
dimensional crystal. Bulk correlation functions probe the 
fluctuations of the bulk degrees of freedom at finite (al- 
beit small) wavevector. The elastic degrees of freedom 
of the lattice or supersolid are the two components of 
the lattice displacement (phonons), u x (q) and u y (q) (or 
equivalently u l(c\) and ur(q )). Hence in the definition 
of S ss (q), Eqs. |3.10 and 3.30| , the average is over the two 
lattice displacements (as well as over the defects), not 
the components of the strain tensor, w a i, which can be 
written as derivatives of the displacements and there- 
fore are not independent degrees of freedom. In contrast, 
the response functions, that most directly relate to ex- 
perimental probes of elasticity, measure the response to 
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a (often uniform) stress applied at the boundary of the 
solid. The corresponding q = 0, uniform deformations 
are described by three zero-mode independent degrees of 
freedom, corresponding to three macroscopic strains u xx , 
u yy and u xy that can be independently induced in a solid. 
As a consequence, for a supersolid (or a crystal in gen- 
eral) the equal-time correlation functions differ from the 
corresponding static response functions, as we now ex- 
plicitly demonstrate. 

Finally, we note that the situation is different in a liq- 
uid, where dislocations have proliferated and act as addi- 
tional degrees of freedom. The proliferated dislocations 
lead to multi-valued lattice displacements and thereby 
allow for a transverse part of the strain Wy, which (for 
the first index i), is forbidden in the supersolid. As a re- 
sult, in a two-dimensional liquid, all three components of 
the symmetric strain (u xx , u yy and u xy ) are independent 
degrees of freedom. 

Returning to the derivation of the linear response to 
a perturbation described by Eq 



3.37 



we treat Uij, w 



and the defect densities as independent degrees of free- 
dom and minimize the total free energy, F tot = F ss + F a , 
with respect to them, to obtain 

— - 2c66Uij + 5ij(cil ~ 2c66)u/cfc + Sij—Snd - <7ij , 



no 



= 



SF tot i 

C44,W z j H tdj , 



Sw z 



"o 



= 



SF tot x~\ 7 , 1 x 
ond Uq n 2no 

= 0, 

mot v . , y 

— — = Ktdi-\ w zi , 

ot d i no 

= . 



(3.38) 
(3.39) 
(3.40) 
(3.41) 



Given that the app lied stress a a i s purely within the 
xy-plane equations ( 3.39 ) and ( 3.41 ) give 



tdi 
w zi 



0, 
0, 



(3.42) 
(3.43) 



and the response is z independent. Solving the remaining 
equations for the strain and the defect density, we find 



Ukk 

Sn d 
no 



1 + 7X 



2(cn - c 66 - 7 2 x) 
X(cn - c 66 +7) 



2(cn ~ cm - l 2 x) 
which, when used inside Eq. 3.38| give 

Uij = R s ilkl a kl , 



(3.44) 
(3.45) 

(3.46) 



where i£f? kl is the uniform static response function given 
by 



R 



(2c 66 - cu)(l + 7x) + 27x(cn - c 66 + 7) 



4c 66 (cn - c 66 - 7 2 x) 



S% j Ski 



2c, 



-SikS 



Jl 



(3.47) 



As a check we observe that for \ — > 0, which freezes 
out the defects, Rfj y ki reduces tolhe well-known response 
function for a defect-free crystalcJ 



R 



crystal 
ij,kl 



(2c 66 - en) 



4c 66 (cn 
SikSji ■ 



C66) 



SijSki 



2c, 



(id 



(3.48) 



Consistent with our discussion above, Rfjki differs from 
the corresponding q z — equal-time correlation function, 
computed with the lattice displacements as the indepen- 
dent degrees of freedom. 



Using Eqs.3.46 and 3.47 we can now compute the re- 
sponse to any uniform stress, <7y . For example, a uniform 
hydrostatic pressure Sp corresponds to 



-SpSi 



(3.49) 



and we obtain the bulk modulus for the vortex supersolid 
as 

1 ISA 
~bTs = ~A~5p~ ' 
1 Sn 

no Sp ' 
_ _Ukk_ 1 5n d 
Sp no Sp ' 
B + X - 1 + 2 1 



Bx -1 _ 7 2 



(3.50) 



where B = c\\ — cqq is the standard definition for a 
defect- free crystal bulk modulus. We note the similar- 
ity in form with the correlatio n func tion S' ss (qj^, q z =0), 



that we used to define cfi, Eq.|3.20|. As discussed above, 



the computation of the uniform static response function 
gives the bulk modulus, c\\ — c^e, while the equal-time 
correlation function (at q z — 0) defines the longitudinal 
phonon modulus, C\\. 



As we observed above for the effective c{\ in Eq. 3.24 



here too in the absence of interaction between the defects 
and the lattice (7 = 0), the supersolid bulk modulus is 
simply determined by the bulk moduli of the two systems, 
added in parallel, -B" 1 = B^ 1 + \- 

Signific ant c aution must be applied in comparing ex- 
pression (3.5C) with experiments. If, for example, the 
compression is performed on a time scale that is slow 
compared to the relaxation time of the elastic degrees of 
freedom, but fast compared to that of the defects (which, 
being a conserved "charge" density, relaxes only diffu- 
sively) , then the defects are effectively frozen on the time 
scale of the experiment. In this case, corresponding to 
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X = 0, only elastic degrees of freedom can respond and 
the result for B ss is simply B of the defect-free crystal. 

The vortex supersolid shear modulus c|| can be probed 
by applying a constant force normal to two opposite lat- 
eral boundaries of a rectangular rry-crossection sample. 
The corresponding stress is given by 



(3.51) 



Such stress induces both the longitudinal (ul) and the 
transverse (mt) deformations of the lattice. The shear 
modulus is defined in terms of the normal strain, u xx — 



1 2 ^"yy ^xx 



(3.52) 



Consistent with our earlier calculation that used equal- 
time correlation function, we find here that 



c 6 6 = c 66 



(3.53) 



i.e., the shear modulus remains unrcnormalized by fluc- 
tuations in density of vacancies and interstitials. 



IV. DISLOCATION LOOPS IN A VORTEX 
SUPERSOLID: A MODEL OF VORTEX LIQUID 

We now turn to the description of a vortex liquid. This 
is a disordered, dissipative state of the flux-line array, 
that results from either a direct melting of a defect-free 
vortex crystal or (possibly) a continuous meltingEij of a 
vortex supersolid, discussed in previous sections. 

Following ideas of Kosterlitz and Thouless,L3 extended 
in Ref. [16] to three-dimensional vortex systems, we de- 
scribe a flux-line liquid as a supersolid with a finite 
equilibrium concentration of unbound dislocation loops 
and vacancy/interstitial defect lines. Such an approach 
nicely complements the more conventional hydrodynamic 
description of the vortex liquid studied in Refs. [37 36 
While somewhat more involved, the advantage of the 
approach taken here is that it provides a more direct 
connection between the vortex liquid and vortex ordered 
phases (crystals), in which the defects are bound, thereby 
presenting a unified description. It also provides a valu- 
able detailed "microscopic" characterization of the dis- 
tinction between vortex liquids and solids. 

The properties of a flux-line lattice in the presence of 
an unbound gas of dislocation loops-but no vacancies 
nor interstitials, were studied in Ref. [l(| In the absence 
of vacancy and interstitial defects, dislocation loops are 
constrained to lie in a plane spanned by their Burger's 
vector and the z axis. These planar dislocation loops 
can only relax applied stresses by "gliding" along the z 
axis. Clearly, once dislocation loops proliferate on all 
scales, vacancies and interstitials will also unbind and 
both of these defects will exist in the resulting bond- 
orientationally ordered liquid. The goal of this section 



is to incorporate vacancies and interstitials into a com- 
plete description of a vortex liquid. We will explicitly 
demonstrate that a finite concentration of vacancy and 
interstitial defects allows for "climb-like" distortions of 
dislocation loops, which can move out of the z-b plane 
by absorbing and emitting vacancies and interstitials. 



A. Model 

To construct a complete model of a vortex liquid we 
now proceed to incorporate dislocation loops into the 
model of a supersolid, presented in Sec. III. We do this by 
allowing multi- valued lattice displacement fields, u. The 
loop integral of u, 



dui(r) = -bi(r) 



(4.1) 



enclosing a dislocation line, fails to close by a Burger's 
lattice vector, b. The direction of integration around the 
contour is that of a right-handed screw advancing paral- 
lel to a unit tangent vector r to the dislocation line. The 
peculiarity of dislocation in a lattice of z-directed lines 
is that while the Burger's vector is two-dimensional and 
by definition lies in the xy plane, the tangent r to the 
defect line is a three-dimensional vector .EjO 

To study properties of the system on scales that are 
long compared to the spacing between dislocation lines, 
we use a continuum description. We consider a small hy- 
drodynamic volume and introduce the Burger's "charge" 
density tensor apj(r), whose integral over an open sur- 
face S, gives the total Burger's vector of dislocation lines 
directed along the surface normal np and enclosed by a 
contour C bounding the surface, 



" s 



(4.2) 



For a single dislocation line, directed along the tangent 
Tp, with Burger's vector the defect density tensor is 
given by ctpi — rpbid^ 2 ' (r±). The rectangular (3x2) den- 
sity tensor apj(y) is therefore a measure of the number of 
dislocation lines with Burger's vector bj crossing a unit 
area normal to the dislocation tangent Tp. We remind 
the reader that Roman letters i,j,k,... are used to de- 
note indices that run only over the values x and y, and 
Greek letters a,/3, 7, ... are reserved for indices that run 
over the three-dimensional set x, y, z. 

By definition, dislocations in the z-directed line crystal 
must have their Burger's vectors lie in the xy plane. Con- 
sequently the three-dimensional vector a Z j(r), (J = x,y) 
describes z-directed edge dislocations. Edge dislocations 
lying in the xy plane are described by the antisymmetric 
part of the two-dimensional tensor CKy(r), (i — x,y and 
j = x,y, with apj = (a^,a Z j)). In the absence of va- 
cancies and interstitials, however, these type of xy-plane- 
directed edge dislocations correspond to a branching or 



12 



merging of flux lines, which necessarily involve fractional 
or double flux- "charged" vortex lines, both energetically 
forbidden. In contrast, vacancies and interstitials allow 
for this type of edge dislocations, as we will demonstrate 
below. Screw dislocations, on the other hand always run 
normal to flux lines, i.e., they lie in the xy plane, and lead 
to entanglement of the vortex lines .c3 They are described 
by the symmetric part of ai(r). 



By rewriting Eq.4.1 



in differential form and then av- 
eraging the resulting equation over a hydrodynamic vol- 
ume, containing many dislocation lines, we relate the dis- 
location density tensor apj(r) to the local lattice strain 
wpj(r), 



e af}l df3W lk {v) = -a ak (r) 



(4.3) 



Hence in the presence of dislocations (ctpj ^ 0), the strain 
tensor wpj contains a singular part that cannot be writ- 
ten as a gradient (dp) of a single- valued displacement 
Finally, dislocation loops must either close or 



field 

terminate at sample boundaries. This amounts to the 
condition 



dpa 0j (i 







We now proceed to incorporate the dislocation degrees 
of freedom into the model of the supersolid studied in 
Sec. III. We first recall that the fluctuations of the lo- 
cal magnetic induction <5B = (B±,5B Z ) are related to 
the changes in density and orientation of flux lines, with 
contributions from both local strains and vacancy and 
interstitial defects. In the long wavelength limit these 
rela tion s ar e ob tained by inse rtin g 5n and t given by 



Eqs.3.6 and 3.7 in Eqs.2.4 and 2.5 



SB Z = (f> [-noWu(r) + Sn d (r)} , 
SB ±l = (j}o[n w zi (r) + t dl (r)} . 



(4.5) 
(4.6) 



In the presence of dislocations, the strain tensor w a j is 
given by the sum of a regular, longitudinal (on first, 
Greek index) part, defined in Eq.2.3 in terms of the 



derivatives of a single- valued displacement field Uj , and a 
singular, transverse (on the first, Greek index) part, w^i, 
due to dislocations 



= d a Uj 



a-3 



(4.7) 



B y im posing the V • B = condition and using Eqsj4 
and 4.6, we obtain 



d z 5n d 



td = -no(-d z wl 



, ,sing\ 



(4.8) 



It is important to note that, in contrast to the defect- 
free crystal and the supersolid, where the elastic part wpj 
(dpiij) identically satisfies V • B = 0, here this condition 
imposes a nontrivial constraint that couples dislocations 
and vacancy-interstitial defects. This becomes apparent 
by using Eq.[lJ^ to eliminate wpj from the constraint in 
favor of the dislocation density tensor, with the result 



d z Sn d + Vj_ ■ t d 



n e 



aj-Jr) . 



(4.9) 



This important condition is one of the main results 
of our work and is an essential ingredient in the com- 
plete "elastic" description of the vortex liquid state. 
We first note that in the supersolid phase, where dis- 
location loops are bound and lattice displacements are 
singe-valued, apj = 0, and the constraint reduces to 



the vacancy-interstitial line continuity condition, Eq.3.8 



On the other hand, in a description of a vortex liquid 
which ignores vacancies and interstitials it reduces to the 
e ij a i] = condition of Ref. 16, enforcing the constraint, 
discussed above, that in the absence of vacancies and 
interstitials, dislocation loops are confined to lie in the 
plane defined by their Burger's vector and the average 
external magnetic field. In the presence of vacancies and 
interstitials Eq.4.9 explicitly demonstrates that these de- 



fects provide a mechanism by which the dislocation loop 
can effectively climb out of this plane. They do this by 
emitting or absorbing vacancy and interstitial line de- 
fects. The allowance of nonplanar dislocation loops in 
vortex supersolids and liquids then removes the artifi- 
cial condition of a vanishing antisymmetric part of the 
(4.4) dislocation density tensor and in general tijCtjj ^ in 
these phases. As anticipated above, this allows xy-plane- 
directed edge dislocations to pro liferate upon melting. 

The above constraint, Eq.4.9, is one way that the ad- 
ditional degrees of freedom associated with proliferated 
dislocations enter the description of the vortex liquid. 
Dislocations of course also contribute directly through 
the "elastic" free energy in a way that we now derive. 

The free energy F of a lattice with dislocations and 
vacancy /interstitial defects is given by F ss = Fi at t + F d + 
Fi nt , but with derivatives of the pho non field replaced by 
the total strain tensor given in Eq. 4.7, i.e., 



+ i / dr [x^(^d) 2 +K(t d ) 2 ] 



(4.10) 



H / drHwudrid + ^'w zi t di \ 

n J 

To make the notation more compact, the purely elastic 
part of the free energ y, giv en by the first term on the 
right hand side of Eq. 4.10, has been written in terms of 
the elastic tensor, 

CaiPj — CQd(S a p — S az Sp z )Sij + CQ()S a jSpi 

+ (cn - 2c 66 )S ai Spj + C^SazSpzSij . (4.11) 

We have also added to the free energy of the defective 
lattice a term F core representing the core energy of the 
dislocations, given by 

^ orc = / dr[E e al + E s *i + E>a^ 3 + E^f]. 

(4.12) 
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The core energy F COIe has been written on the basis of 
symmetry considerations. It incorporates terms account- 
ing for the edge and screw dislocation core energies per 
unit length, E e , E' e , E s and E' s . As discussed in Ref. 
fL6|, although E' e = E' s = for a single dislocation line, 
nonzero values of E' e and E' s are required to describe short 
range interactions in the hydrodynamic limit. The val- 
ues for the core energies are estimated to be_E e ~ c§§b 2 , 
E s ~ E' s ~ ^mb 2 and E' e ~ 7cTT^I& 2 .N 

While the elasti c pa rt (d a Ui) of the strain tensor iden- 
tically satisfies Eq.4.3, the singular part is found by "in- 
verting" this equation in Fourier space. This gives 



(q) 



e a p 7 qpa J3 {q) +iq a ipj{q) , (4.13) 



where tpj is an arbitrary function, reflecting the fact 



that the solution to Eq.4.3 is only determined up to 
an arbitrary longitudinal part. This is analogous to 
the gauge freedom that appears in electromagnetism (or 
other gauge theories) when, for example, Maxwell equa- 
tion V • E = p is solved for the electric field E in terms 
of the charge density p. 

One convenient and natural choice of ipj is obtained 
by requiring that the corresponding xif a ^ minimizes the 
total free energy F, Eq. 4.10 , with respect to lattice dis- 
placements. Such a choice is mathematically convenient 
for computing thermodynamic averages, because it re- 
moves all couplings between the phonons Uj and the de- 
fect degrees of freedom. Physically, this choice of to^ • 
corresponds to elastically equilibrated defects. 

To compute ipi correspo ndin g to this con venie nt gauge 
choice, we insert w a i, Eq.4.7 into F, Eq.4.10, and re- 



quire that the linear terms in u vanish, which of course 
is equivalent to the requirement that uf ai minimizes F. 
The resulting Euler-Lagrange equation is given by 



SF 
5u 



7 i 
— di6n d d z t d i 



0, 



"o 



n 



where 



= C ai 0jWf3j 



(4.14) 



(4.15) 



is the corresponding stress tensor. The solution is more 
conveniently written in Fourier space, where it is given 
by 



F 



dof 



n 

1 

2 



1 ) ij la C aJ f3k e/3 7 r;97 a r)k (q) 

dq 



(4.16) 



(2tt) 3 



Here A 1 is the inverse of a 2 x 2 matrix A, with 
iPjQ/3- Its elements are given by, 



1 



r T (q) 



(gu ~ c 6 6 ) 
Ti(q) 



-q±iq± 



(4.17) 



where IV (q) and r^(q) are the transverse a nd lo n gitu- 



dinal elastic elastic stiffnesses, defined in Eqs. 3.16, 3.17 



With this choice of ipi, by construction, the total free 
energy breaks up into two parts 



F = F 



latt 



F, 



del, 



(4.18) 



where Fi at t is the defect-free elastic part of the free en- 



ergy given in Eq. 2.1 and rewritten here for convenience 
in Fourier space, 



F l!x 



d 3 q 

(2tt) 3 



r T (q)Mq)| 2 +T L (q)\u L (q) 



(4.19) 



This involves only phonons degrees of freedom, as for a 
defect-free crystal. The second, defect part, Fdef, is the 
free energy of an interacting gas of dislocation loops and 
vacancy-interstitial line-defects. It is given by 



-Fdef 



drw s al C t 



2n 2 
1 

n 



dv[x~\5n d ) 2 +K{t d ) 2 } 
dr[yw° i 8n d + j'w s zi t dl ], 



(4.20) 



where 
by Eqs. 



,- is t he si ngula r part of the strain tensor, given 
4.13| and |4 . 1 6|. This defect free energy must 

4~£ 



be sup plem ented with the continuity conditions Eq 
and Eq.4.4 for the vacancy/interstitial defect lines and 
the dislocation lines, respectively. The constraint (4.9) 
can be directly incorporated into the free energy by us- 
ing it to explicitly eliminate the longitudinal part of the 
vacancy-interstitial defects tangent vector, 



t 



d(q) = q-L •trf(q) » 

Qz 



q± 



-5n d (q) 



-no£ya,,(q) 



(4.21) 
(4.22) 



in favor of the defect density Sn d (q) and the antisymmet- 
ric part of the dislocation density tensor (q) . Finally, 
using Eqs. 4.13 and 4.16 to eliminate the singular s train 
field in terms of the defect degrees of freedom in Eq 
we obtain 



4.2C 



^wi(q)a^(q) a fi(-q) + A(q)5n d (q)5n d (-q) + C(q)t d (q)f d (-q) 



+ iD f _ lt (q) a fli (-q)5n d (q) - a Ali (q)5n d (-q) +iG^(q) a^(-q)tj(q) ■ 
where tj(q) = t d (q) — q±t^(q) is the transverse part of the defect tangent vector. 



i(q)*d(-q) 



(4.23) 
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The kernels ^^-(q), A(q), C(q), D^(q) and G Mi (q) 
depend in a complicated way on the elastic moduli 
(c66, en, C44) of the lattice, and on wavevectors q± and 
q z , and are given explicitly in Appendix B. 



B. Correlations and response functions 

In this section we evaluate the renormalization of the 
elastic constants of the flux-line lattice due to dislocations 
and vacancy/interstitial defects. T he rcnorm alized elas- 
tic constants are defined by Eqs. (3.20-3.22), but with 
the understanding that the correlation functions are now 
those of a lattice with an equilibrated concentration of 
defects and dislocations. This means that the structure 
function and the tilt correlation function are formally 
given in terms of the str ain te nsor and the defect fields 
by the same expressions 3.10| and 3.11 used for the su- 
persolid, but the strain tensor w a i is now the total strain 
given in Eq. 4.7, including the singular part. The same 



holds for the correlation function that determined the 
shear modulus defined on the first line of Eq. p. 22 . The 



brackets in these correlation functions now denote a ther- 
mal average with t he free energy Fi at t + -Fkof , with fdof 
given by Eq. 



4.23 



The average is carried out by inte- 
grating over all configurations of dislocations, described 
by the components of the dislocation density tensor, api, 
and vacancy/interstitial lines, described by the defects 
density, Srid, and tilt field, t^. The integration must be 
done subject to the constraint that dislocation lines are 
closed, given by Eq. 4.4. The "continuity" constraint 
for defect lines, expressed by Eq. 4.E, has already been 
incorporated into the free energy i<def ■ The computation 
of correlation functions is conceptually simple but tech- 
nically quite involved and was carried out using Mathe- 
matica symbolic manipulator program. Only the results 
will be given here. 

The full expression for the density, tilt field, and other 
correlation functions are too "horrifying" to be shown 
here, and we therefore only display their long wave- 
length limits, which determi ne the rcn ormalized elastic 
constants, according to Eqs. ( 3.20 - 3~22| ). 



The structure factor vanishes as q± — ► 0, as required 
by the density sum rule. The renormalized longitudinal 
modulus as defined by Eq. ( p. 20 ) is given by, 



1 

Cll 



cii[2c 66 (cn - c 66 ) + CnE e q\] + 2c| 6 x 1 + 7(7 + 2cn)(2c 66 + E e q\) 
X _1 [ 2c 66(cn - cm) + cnE e q\] - 7 2 (2c 66 + E e q\) 



(4.24) 



The renormalized tilt modulus defined by Eq. (3.21) 



is 



1 ™ / c 44 + K - 2-f' 

= hmT T (q z , qi _ = 0) = „ J . (4.25) 



A) 



C44K — 7 /2 



We find that in the long wavelength limit both and 
C44 are identical to the elastic constants of a lattice with 



only vacancy and interstitial defects. In other words, 
somewhat surprisingly, the coupling of dislocations to 
vacancy /interstitial defect lines does not yield any addi- 
tional renormalization of the tilt and compressional mod- 
uli, even if it does make it possible for dislocations to 
relax stresses by climbing out of the (b, z) plane and al- 
lows for (otherwise forbiddingly costly) xy directed edge 
dislocations. Finally, the renormalized shear modulus is 



1 



1 



X 1 cii(cn - 2c 66 ) - 7 2 (cn - 4c 66 ) 



(id 



(ffx) 



C66 



E eV± en X _1 ( 2c 66(cn - c m ) + cnE e q\) - 7 2 (2c 66 + E e q\ 



(4.26) 



Dislocations renormalize the long wavelength (q — > 0) 
shear modulus to zero, yielding liquid-like response to 
shear stresses. 

We now discuss various limiting cases for our re- 
sults. In the absence of coupling between the gas of va- 
cancy/interstitial defects and the lattice (7 = 0, 7' = 0), 
the various correlation functions are simply given by the 
sum of the contributions from a lattice with an equilib- 
rium concentration of unbound dislocations (correspond- 
ing to the constrained, hexatic line liquid discussed by 
Marchetti and NelsonES) and from a liquid of defect lines. 
The corresponding elastic constants add in parallel, with 



(cfi(<7i-)) 



(<?_!_)'' 7=7'=0 



x + 



JIN • 

--11 



(4.27) 



-44(31 
1 



) / 7 = 7 '=o 



1 

K + 
1 

JVIN : 



JV1N ■ 

-44 



(4.28) 
(4.29) 



-66Vy-U ' 7=7' =u i- 66 

where we have denoted by the superscript "MN" the elas- 
tic constants of a lattice with an equilibrium concentra- 
tion of unbound dislocations, but with ^constraint for- 
bidding vacancy and interstitial defects, £3 given by 



„MN 



11 



G 66 



C44 



1 

Cll 

1 

C66 
1 

C44 



1 



2c 2 

zc 66 



2c 6 6(cn 

1 



C66) + cnE e q 2 j_ 
(en - 2c 6 6) 



E e q\ 2c 66 (cii - c 66 ) + c n E e q2_ 



(4.30) 



(4.31) 
(4.32) 
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The elastic constants of a supersolid, i.e., the lattice with 
only vacancy and interstitial defects, can be obtained 
by letting all the dislocation core energies go to infin- 
ity (E e , E' e , E s , E' s — > oo). It is easy to see that in this 
limit we recover the results discussed in Sec. II. Finally, 



when x 



and K — + oo vacancy and interstitial de- 



fects are forbidden and the clasti c constants reduce again 
to those given in Eqs. p3o[]4~32l . 



We recall that the coupling of dislocations to vacancy 
and interstitial lines allows dislocations to relax stresses 
by climbing out of the (b, z) plane via the emission or ab- 
sorption of vacancies and interstitials - a mechanism for 
relaxing stresses that isi forbidden in the hexatic liquid of 
Marchetti and NelsorJla. 



* 
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Fig. 5 . Gliding of an edge dislocation loop, allowed even in 
the absence of vacancies and interstitials. For simplicity 
we have sketched the dislocation loop for the case of 
a square lattice. The dotted vertical lines refer to the 
defect-free square lattice and serve as a guide for the eye. 



Figure 5 shows a planar dislocation loop withJp = 6x 
of the type considered by Marchetti and Nelsor£3. This 
loop lies in its glide plane (the xz plane) and can easily 
relax a shear by gliding in this plane. The climbing of the 
same dislocation loop out of its glide plane is described 
pictorially in Fig. 6. 




Fig. 6. Climbing of an edge dislocation loop in a square 
lattice, only allowed in the presence of vacancies and 
interstitials. The Burger's vector of the loop is in the +x 
direction. The lattice configuration at two cross-sections 
z\ and Z2 shows the climbing in the y direction, normal 
to the Burger's vector. 

It is clear that climb can only occur via the emission 
or absorption of vacancy or interstitial defects, as sum- 
marized by Eq.I.E. Such climb "motion" can occur in 



response to the force on the dislocation loop resulting 
from applying a uniform external tilt to the lattice, cor- 
responding to a constant stress of** = Si x 8 az a^. The 
force on a dislocation line due to a uniform stress of** is 
the familiar Peach-Kohler force, given by 



pPK _ ,-rCXtl 

r f3 — t af3jT^a ia l>i 



(4.33) 



Equation 4.33 differs slightly from the corresponding ex- 
pression found in textbooks as the stress tensor in a flux- 
line lattice is not a symmetric (or even a square) matrix. 
The Peach-Kohler force on the rectangular loop shown 
in Fig. 5 due to a uniform tilt in the xy plane is normal 
to the plane of the loop (along y) and there is no force 
on any segment of the loop parallel to z. Specifically, for 



= x, we find F 



PK.tilt 



by and for r 



-x, we 

find F PK,tllt = af^by. In other words, the Peach-Kohler 
force acts as a couple and tries to rotate the loop out of 
its glide plane. The screw components of the loop run- 
ning along ±x can glide in any plane that contains them. 
In particular, they will glide in the y or — y direction un- 
der the action of the force. For this to happen the edge 
sections must climb out of the glide plane by emitting or 
absorbing vacancies and interstitials, as shown in Fig. 6. 
By gliding out of the sample along the y direction, the 
screw dislocations can relax a uniform tilt of the lines 
towards ±x. We note that if o ne thinks of z as a fic- 
titious time, the constraint 4.9 is formally identical to 



the temporal continuity equation for the density of point 
vacancy /interstitial defects, diffusing inr-the presence of 
dislocations in a two-dimensional latticeEa Dislocations 



1G 



provide a source of point defects as they climb across 
the sample. The renormalization of the tilt modulus can 
therefore occur only when vacancy and interstitial defects 
are allowed. 



It is also instructive to consider the behavior of another 
typical dislocation loop, shown in Fig. 7. This loop has 
b = fox and lies in the yz plane. The two edge segments 
running parallel to z lie in two different and parallel glide 
planes. The loop is closed by segments running along y 
which are also edge dislocations in nature, as b _L t. 



This is an example of a dislocation loop characterized 
by Cijajj 7^ 0. Such a loop is not allowed in the absence 
of vacancy and interstitials, as without such defects the 
edge segments running along y would require a row of 
vortex lines to merge or split into single lines carrying 
twice or half a flux quantum - an energetically forbidding 
configuration. Denoting by z\ < Z2 the vertical location 
of the two segments of the loop running along y, we can 
understand the existence of the loop as arising from a set 
of interstitial defects that are randomly distributed in the 
xy plane for z < Z\ and z > z 2 > but organize themselves 
into a vortex sheet for z\ < z < Z2, acting like an extra 
row of lines in this region. In the region near z\ and z%, 
clearly d z nd ^ 0, corresponding to a nonvanishing value 
of 




_L 



Fig. 7. An edge dislocation loop described by a nonvan- 
ishing afj6ij, allowed only in the presence of vacan- 
cies and interstitials. The bottom figure showing a side 
view of the loop emphasizes the equivalence between va- 
Under a constant stress of? from a uniform tilt cancy/interstitial defects and a nonvanishing a^ey. The 



applied to the system, the Peach-Kohler force on the loop 
consists again of two forces of equal magnitude applied 
on the sections running along y and directed along ±x, 
as there is no force on the sections of the loop running 
along z. Under the action of this couple, the loop rotates 
out of its plane. While this motion occurs in the glide 
plane for each section of the loop, it requires "diff usion" 
of vacancies and interstitials according to Eq. 4.9. 



loop can be thought of as arising from an insertion of 
a row of "finite- length vortices", which correspond to a 
wandering of interstitial defects. 

The addition of dislocations also removes the nonana- 
lyticity of the transverse part of the tilt-tilt correlation 
function, present in both the defect-free vortex lattice 
and the vortex supersolid. In fact we obtain 



Tr(q±,qz = 0) = n\k B T 



K + C44 + 2(E S + E' s )q\ — 2j' 
KC44 + 2K(E S + E' s )q\ - 7 ' 2 ' 



(4.34) 



and 



(W 2 



lim Tt(qa 

n\k B T 

c 44 



, q* = 0) 



(4.35) 



with C44 given by Eq. i.2l. As indicated in Eq. [2.23 
lim^^o Tr(q±, q z = 0) determines the winding number 
(W 2 ) and the corresponding superfluid density n s of the 
equivalent boson system, (W 2 ) = %n s /m. The analytic- 
ity of the transverse tilt-tilt correlator at long wavelength 
indicates that n s = no, i.e., the bosons are in the super- 



fluid state. Conversely, this corresponds to an entangled 
flux-line array, with no longitudinal superconductivity. 

For comparison, in the model of a constrcmed hexatic 
line liquid, studied by Marchetti and Nelson,li3 where the 
absence of vacancy /interstitial defects prevents disloca- 
tion loops from climbing out of the (b, z) plane, one ob- 
tains 



T^(q ± ,q z = 0) 



Cii + 2{E S + E' s )ql 



(4.36) 



and 
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lim T^ N (gj 



,<lz 



0) 



n\k B T 



c 44 



(4.37) 



We find that independent of whether or not vacancy and 
interstitial defects are included in the description of the 
hexatic vortex liquid, the transverse tilt-tilt correlator is 
analytic and lim gj ^o 2t(<?_L; <Zz = 0) is finite, indicating 
that both systems are entangled. Vacancies and intersti- 
tials do, however, decrease the tilt modulus, further en- 
hancing the entanglement of the vortex array. Although 
vacancy and interstitial defects must be incorporated in 
a consistent description of a flux-line hexatic, they do 
not change the qualitative properties of the hexatic liq- 
uid, which even in their presence maintains a vanishing 
long- wavelength shear modulus and does not exhibit lon- 
gitudinal superconductivity. 

To further characterize orientationally ordered liquid 
it is useful to define a characteristic length scale £j_ , that 
determines the typical transverse size of a disentangled 
flux-line bundle. This is a region of the flux array where 
the lines remain locally disentangled in the limit of infi- 
nite sample thickness along the field (z) direc tion. In the 
absence of vacancies and interstitials, Eq. ( 4.36 ) gives 
~ ^/2(E S +E' s )/ Cii . Using E a ~E',~ ^E^b 2 , 

we find ~ ao^cee/cii} 1 ^, with clq = ^/l/no the 

mean intervortex separation. In flux-line arrays we typi- 
cally have C66 << C44 and the flux-line hexatic is entan- 
gled over all macroscopic scales. 

The coupling of dislocations to vacancy and interstitial 
defects renor maliz es this entangling length £j_ , as can be 
seen from Eq.4.34, which for small q± is given by 



r T (?x,o) 



nlk B T 



2 c2 



with 



-44 



\K 



2{E S +E' S ) \c ii +K-2 1 '\ 



(4.38) 



(4.39) 



It is easy to see that for all physical parameter values, 
£-L < £1™ ■ Not surprisingly, we find that vacancies and 
interstitials therefore decrease the typical size of a disen- 
tangled flux bundle. 

Finally, we have also studied the effects of vacancy and 
interstitial defects on the properties of the orientationally 
ordered hexatic liquid. Although the finite wavevector 
behavior is considerably modified, we find that at long 
wavelength the effective hexatic stiffness is still given 
by the expression obtained in Rcf. 16. This result is 



consistent with the lack of long-scale orientational order 
and therefore a vanishing orientational stiffness in the 
isotropic line-liquid of vacancies and interstitials. 



V. CONCLUSIONS 

In this paper we have studied the effects of vacancy 
and interstitial line defects in the supersolid and the liq- 



uid phases of flux-line arrays in the mixed state of type-II 
superconductors. 

The transition from a fully ordered Abrikosov crystal 
phase into the vortex supersolid phase takes place at a 
critical temperature or magnetic flux density, at which 
vacancy and interstitial defects proliferate, providing a 
mechanism for vortex entanglement and yielding finite 
resistivity. In order to study the long wavelength elastic 
properties of the supersolid phase, we modeled the super- 
solid as a lattice with an equilibrium concentration of un- 
bound vacancy and interstitial defect lines and computed 
the renormalization of the elastic constants due fluctu- 
ations in the density and orientations of these defects. 
As expected, at long wavelengths, we find that vacancy 
and interstitial density fluctuations do not renormalize 
the shear modulus, confirming the finiteness of the su- 
persolid shear modulus, which distinguishes it from the 
vortex liquid. In contrast, these defect fluctuations yield 
a finite downward renormalization of both the compres- 
sional and tilt moduli. The renormalization of the tilt 
modulus stems from the fact that the liquid of vacancy 
and interstitial defects promotes flux-line wandering. We 
explicitly demonstrate that this defect-based vortex line 
derealization mechanism spoils the Meissner response to 
a shear tilt perturbation, characteristic of the defect-free 
vortex lattice. In other words, the vortex supersolid does 
not exhibit longitudinal superconductivity and, in this 
respect, similarly to the vortex liquid, is always entan- 
gled. 

It is clear that the vortex liquid phase, where unbound 
dislocation loops have proliferated, is also characterized 
by a finite density of vacancy and interstitials, which 
must, therefore, be included for a consistent description 
of such a resistive phase. A model of the vortex- line liquid 
as a solid with an equilibrium concentration of unbound 
dislocation loops, but no vacancies and intcpstitials, was 
previously studied by Marchetti and NelsonllS. As a con- 
sequence of no-vacancy/interstitial constraint implicit in 
their model, only planar loops, that cannot climb out of 
their glide plane, were included. We have generalized this 
model by allowing for a gas of vacancy and interstitial line 
defects, coupled to the dislocation loop gas and the elastic 
(phonon) degrees of freedom, to be present in the orien- 
tationally ordered vortex liquid. We explicitly demon- 
strated that a finite concentration of these vacancy and 
interstitial defects allows the dislocation loops to climb 
out of their glide plane by emitting and absorbing vacan- 
cies and interstitials, and allows for new types of edge 
dislocations (otherwise forbidden), thereby significantly 
increasing the entropy of topological defects characteriz- 
ing the vortex liquid. T his physical result is mathemati- 
cally summarized by Eq. |4.S| , in which non-planar disloca- 
tion loops and horizontal edge dislocations, with Burger's 
vector in the xy-plane, act as sources for vacancy and in- 
terstitial line defects. This resulting climb-like motion 
provides a mechanism for the relaxation of an externally 
applied tilt, not present in the model of Ref. nq. We 
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find, however, that at long wavelengths there is no in- 
dependent renormalization of the tilt and compressional 
moduli by dislocations, other than that already induced 
in the supersolid by fluctuations in the density and ori- 
entation of vacancy and interstitial defects. This result 
is somewhat puzzling and may be a consequence of the 
Debye-Huckel approximation used to here. We are cur- 
rently investigating the extension of the model presented 
here beyond the (quadratic) Debye-Huckel treatment of 
defects in order to answer this question. The coupling 
of dislocations to vacancy and interstitial defects does, 
however, affect the response to a shear tilt perturbation 
that probes longitudinal superconductivity. While both, 
the vortex supersolid and the vortex liquid fail to exhibit 
longitudinal superconductivity, and are therefore always 
entangled, in the supersolid it is only the response of the 
vacancy and interstitial defects (which form only a small 
fraction of the total vortex-flux density) that spoils the 
Meissner effect, as 



lim Xt (l±,Qz 



0) 



4vr V 



D 



inKJ 



X 



(5.1) 



where K is the tilt modulus of the vacancy-interstitial de- 
fect gas. The transverse susceptibility \t i while larger 
than — 1/47T- the value required for a perfect Meissner re- 
sponse - remains negative (diamagnetic) and finite. We 
expect therefore still an appreciable screening of longitu- 
dinal currents, running along the applied field direction. 
In contrast, for the orientationally ordered hexatic vortex 
liquid (where dislocations loops also unbind), we find 



lim n x¥ x (q±,Qz = °) = -7-( 1 - 

q±^0 4-7T V 



D 2 



47T C f 4 



(5.2) 



where is the signi ficant ly reduced vortex liquid tilt 
modulus, given by Eq ]4.25| . At all, but very low vortex 
densities C44 ss B 2 /Air and the diamagnetic susceptibil- 
ity is therefore vanishingly small. Hence, as expected, 
the fully entangled vortex liquid has essentially normal 
response to a shear tilt perturbation and exhibits no lon- 
gitudinal superconductivity. 

We hope that the results for the properties of the su- 
persolid and hexatic phases, and as well as for the tran- 
sitions scenarios from and into them, presented here, 
will serve as a usefulpguide for identifying these exotic 
phases in experimentsE-3 and simulations^ of vortex sys- 
tems. It is also likely that the equilibrium analytical 
description developed here will be useful for incorpo- 
rating defects . into , the current elastic theories of driven 
vortex latticescacj, as well as for interpreting and char- 
acterizing fhj? J esuits of experimentsc3 and numerical 
simulations 4 Ir^ on these rich nonequilibrium systems. 
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VII. APPENDIX A: DETAILS OF COUPLING 
CONSTANTS APPEARING IN THE 
SUPERSOLID MODEL 

In this appendix we study the phenomenological cou- 
plings that appear in our model of a sup ersolid, defined 
by the effective free energy F as , Eq. |3.12 . 



Our goal is to 

determine the range of values that these parameters can 
ass ume in a physically realistic model of a supersolid. In 
Sec. Ill B we found that the effective supersolid moduli 
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and C44, Eqs. 3.23,3.24 vanish at intermediate values 



of the couplings 7 and 7', suggesting a vortex density 
instability within the supersolid phase. The purpose of 
this appendix is to show that such an instability is only 
apparent and arises from our definition of parameters. 
As we have seen in Sec. [II A, the supersolid effective 



free energy consists of the elastic lattice part, -Fiatt, the 
vacancy-interstitial defect part, F^, and a coupling be- 
tween the elast ic and defect degrees of freedom, Fi n t, 
given in Eq.3.9. The form of Fi nt was dictated by sym- 
metry, with 7 and 7' unknown phenomenological param- 
eters. We argue here, however, that in a more realis- 
tic model of a supersolid, the coupling between elastic 
degrees of freedom and defects arises from fluctuations 
in the magnetic induction, <5B, and therefore should be 
written as an expansion in powers of <5B. To lowest 
(quadratic) order, such energetic contribution, when ex- 
pressed in terms of vortex line density fluctuations 5n 
and t, lead to Fi nt given by 



1 

2^g 



dr 



a(Sn) 2 +a'{tf 



(7.1) 



Inse rting the expressions for 5n and t given in Eqs. 3.6 and 
3.7 in Fi n t, expanding the resulting expression and com- 
bining it with Fiatt and F d , we find F ss = Fi att + F d + F int 



dr 



2c 6 6U 2 . + (en - 2c66)wL + C44(<9 z u) 2 + 



-XjiSnd) 2 



K 



27 

n 



Snd'V ■ u 



2f 
n 



(7.2) 
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Clearly Eq. 7.2 has the same functional form as F ss stud- 



ied in the main text. The coupling constants, however, 
are given by 



Cll 

C66 
C44 

r 1 

k 



en + a , 

C66 , 

C44 + a' , 
X -1 + a , 
K + a' , 

-a , 



(7.3a) 
(7.3b) 
(7.3c) 
(7.3d) 
(7.3e) 
(7.3f) 
(7-3g) 



and they all depend on the two independent parameters 
a and a', which are always positive. 



We can now reexpress the re norm aliz ed su persolid 
moduli c{\ and C44, given by Eqs. 3.23 and 3.24 , respec- 
tively, in terms of the couplings of the model defined here, 
with the result 



-44 



CnX" 1 + (en + X" 1 )^ 

en + AT 1 
c Ai K + (c 44 + K)a' 

C M + K ' 



(7.4a) 
(7.4b) 



which clearly do not vanish (or diverge) for any pos- 
itive values of a and a'. Hence the instability found 



in Sec. HI B was spurious, an artifact of expressing our 



results in terms of 7 and 7' and allowing these coupling 
constants to access values that are unphysical in a generic 
model of a supersolid. 



VIII. APPENDIX B: THE DEFECT FREE ENERGY 



Here we give the expressions for the various kernels contained in the defect free energy of Eq. 4.23. The nonlocal 
kernel R a i,/3j is given by 



Rai,f3j(q) — -^Bai.pjiq) + A a j )( 8j-(q) + E a i£j, 



where 



(8.1) 



2) 



Bai,/3j(<i) = C a i0j — C Q i 7 fc(/ 7 (A 1 )kiqr 1 C rl lf3j e a Aju9A e /3fi/$f j 

describes the long-range interaction between dislocation loops in the absence of vacancy and interstitial defects, while 

(8.3) 



^ai.i/jW — ~^T\ K p J e zai£z/3j ^ — l e zk a q±iezf3j + &zk/3q±j£zai\ 



7 ffz^^g-Lfc 



(8.4) 



is the part of such interaction mediated by vacancy and interstitial defects. The matrix E a i,@j describes the dislocation 
core energy and is given by 

E>aij3j — 2E e 5 az 6i3 Z dij + 2E s S a iSfjj + 2E' s 5 a p(l — 5 az 8p z )5ij + 2E e e za ie z ^j. (8.5) 

The effective interaction between the vacancy and interstitial defects is described by the two scalar kernels 



1 ' K^-^-iqljql^ 
1± 



2 



and 



<?(q) = - 



K 



/2 2 

7 qi 



(8.6) 



(8.7) 



Finally, the tensors D ai and G a i describe the coupling between the dislocation loop gas and the liquid of vacancy- 
interstitial defect lines, and are given by 

Ari(q) = — { [ySpi ~ l'S f 3 Z q± l q z /q±\ - ^- (7 - l'q z /qj_)q fl qxiC t j,ji3i\e^ a q J - — \— + ^-(7 - Yq z /qj_)] e zai , (8. 

Una y- 1 f, > fin I " 1 ' 



G m (q) 



n q 
n q 



l5j3 Z Ci — -^q^CjC^jfii 

1 T 



qi r L 



(8. 



with Ci = ey<?_i_j and q^ = q$/q, 



20 



1 D. Huse and L. Radzihovsky, in Proceedings of 1993 Al- 
tenberg Summer School, Fundamental Problems in Statis- 
tical Mechanics VIII, edited by H. van Beijeren and M. H. 
Ernst (Elsevier, Netherlands). 

2 G. Blatter, M.V. Feigel'man, V.B. Geshkenbein, A.I. 
Larkin, and V.M. Vinokur, Rev. Mod. Phys. 66, 1125 
(1994). 

3 E. H. Brandt, Rep. Prog. Phys. 58, 1465 (1995). 

4 D.R. Nelson, Phys. Rev. Lett. 60, 1973 (1988); D.R. Nelson 
and S. Seung, Phys. Rev. B 39, 9153 (1989). 

5 As first pointed out in Ref. one also expects a vortex 
liquid state to exist just above the H c i(T) curve, in very 
clean superconductor samples. 

6 M. P. A. Fisher, Phys. Rev. Lett. 62, 1415 (1989); D. S. 
Fisher, M. P. A. Fisher, and D. A. Huse, Phys. Rev. B 43, 
130 (1991). 

7 A. Larkin, Sov. Phys. JETP 31, 784 (1970); A.I. Larkin 
and Y.N. Ovchinnikov, Sov. Phys. JETP 38, 854 (1974). 

8 T. Giamarchi and Le Doussal, Phys. Rev. Lett. 72, 1530 
(1994); M.J.P. Gingras and D.A. Huse, Phys. Rev. B 53, 
15183; D. S. Fisher, Phys. Rev. Lett. 78, 1964 (1997). In 
ultra-clean superconductors, the distinction between the 
conventional ordered Abrikosov lattice and the Bragg glasft 
appears, however, only beyond a very large Larkin scale.u 
Analogous topologically ordered, but elastically disordered 
phases also appear in other randomly pinned periodic sys- 
tems, such as for example smectic liquid crystals confined 
in aerogel: L. Radzihovsky and J. Toner, Phys. Rev. Lett. 
79, 4214 (1997); B. Jacobsen, K. Saunders, L. Radzihovsky 
and J. Toner, unpublished. 

9 P.L. Gammel, L.F. Scheemeyer, J.V. Waszczak, and D.J. 
Bishop, Phys. Rev. Lett. 61, 1666 (1988). 

10 R.H. Koch et al., Phys. Rev. Lett. 63, 1511 (1989). 

11 E. Zeldov et. al, Nature 382, 791 (1996). 

12 Although a clean vortex crystal appears to melt via a first 
order transition, as first observed in Ref. [ll], a fluctuation- 
driven continuous melting transition is theoretically possi- 
ble, as demonstrated for model systems in Ref. |l^. 

13 L. Radzihovsky, Phys. Rev. Lett. 74, 4722 (1995); L. Ba- 
lents and L. Radzihovsky, Phys. Rev. Lett. 76, 3416 (1996). 

5. A. Ktitorov, B.N.Shalaev, and L. Jastrabik, Phys. Rev. 
B 49, 15248 (1994). 

14 J. M. Kosterlitz and D.J. Thouless, J. Phys. C6, 1181 
(1973); B. I. Halperin and D. R. Nelson, Phys. Rev. Lett. 
41, 121 (1978); D. R. Nelson and B. I. Halperin, Phys. 
Rev. B 19, 2457 (1979); A. P. Young, Phys. Rev. B 19, 
1855 (1979). 

15 D.S. Fisher, Phys. Rev. B 22, 1190 (1980). 

16 M.C. Marchetti and D.R. Nelson, Phys. Rev. B 41, 1910 
(1990). 

17 D.R. Nelson and J. Toner, Phys. Rev. B 24, 363 (1981). 

18 P. C. Martin, O. Parodi, and P.S. Pershan, Phys. Rev. A 

6, 2401 (1972). 

19 M. P. A. Fisher and D. H. Lee, Phys. Rev. B 39, 2756 
(1989). 



E. Frey, D.R. Nelson and D.S. Fisher, Phys. Rev. B 49, 
9723 (1994). 

L.I. Glazman and A.E. Koshelev, Phys. Rev. B 43, 2835 

(1991) . 

E. Frey and L. Balents, Phys. Rev. B 55, 1050 (1997). 

D. T. Fuchs, E. Zeldov, T. Tamegai, S. Ooi, M. Rappaport, 
and H. Shtrikman, Phys. Rev. Lett. 80, 4971 (1998). 

E. M. Forgan et al, Czech. J. Phys. 46, 1571 (1996). 
Since in type I melting two order parameters (vacancy- 
interstitial density and dislocation density) with unrelated 
symmetry order simultaneously, we expect type I melt- 
ing to be 1st order. In contrast, type II melting can 
be continuous, consisting of 3 intermediate transitions: 
solid-to-supersolid, supersolid-to-hexatic, and hexatic-to- 
(isotropic) liquid (see Ref. [l3| ). 

In a vortex liquid dislocations are, by definition, unbound. 
Furthermore, a pair of dislocations forms a vacancy or an 
interstitial defect. We expect therefore that there will al- 
ways be a finite density of vacancy and interstitial defects 
in the liquid phase. 

The relation to the Lame coefficients fi and A, commonly 
used in the study of elasticity of solids is: C66 = M an d 
en = 2(i + A. 

E.H. Brandt and U. Essman, Phys. Status Solidi B 144, 
13 (1987). 

D.S. Fisher, in Phenomenology and Applications of High- 
Temperature Superconductors, K.S. Bedell et al., eds. 
(Addison- Wesley, 1992), p. 287. 

D.R. Nelson and P. Le Doussal, Phys. Rev. B 42, 10113 
(1990). 

More generally the relationship between the hydrodynamic 
fields and the local induction is nonlocal and can be ob- 
tained by solving the anisotropic London equation, as 
shown for instance in Ref. |30| 

T. Chen and S. Teitel, Phys. Rev. B 55, 15197 (1997). 

D. R. Nelson and V. M. Vinokur, Phys. Rev. Lett. 68, 2398 

(1992) ; Phys. Rev. B. 48, 13060 (1993). 

E. L. Pollock and D.M. Ceperley, Phys. Rev. B36, 8343 
(1987). 

P. Benetatos and M.C. Marchetti, to appear in Phys. Rev. 



B fcond-mat/9808270| ). 

L. Radzihovsky and E. Frey, Phys. Rev. B, 48, 10357 
(1993). 

M.C. Marchetti and D.R. Nelson, Physica C 174, 40 
(1991). 

H.M. Carruzzo and C.C. Yu, Phil. Mag. B 77, 1001 (1998). 
A. Zippelius, B.I. Halperin and D.R. Nelson, Phys. Rev. B 
22, 2514 (1980). 

We note in passing that from this expression we see that the 
stability of a defect-free crystal is determined by the two 
conditions: (i) C66 > and (ii) en > C66 (or equivalently 
in Landau and Lifshitz notation fi > and /i + A > 0), de- 
rived from the requirement of a finite zero mode response 
to a uniform stress. This condition is different and more 
stringent than that imposed by the requirement of a finite 
response of the finite wavevector bulk modes (nr( q) and 
mi (q)), obtained from the propagator of F e i Eq.3.12, which 
gives: (i) ce,a > and (ii) en > 0. 

F.R.N. Nabarro and A.T. Quintanilha, in Dislocations in 
Solids, edited by F.R.N. Nabarro (North-Holland, Amster- 



21 



d ^m, 1980), Vol. 5. 

42 E.H. Brandt, Phys. Rev. B 34, 6514 (1986); Jpn. J. Appl. 
Phys. 26, 1515 (1987). 

43 These renormalized moduli can also be obtained by sim- 
ply integrating out the defect degrees of freedom 8n d and 
t d inside the partition function, and identifying the corre- 
sponding effective elastic moduli in the remaining effective 

elastic free energy functional. 

S. Scheidl and V.M. Vinokur, unplublished (cond 



mat/9702014) 



45 L. Balents, M.C. Marchetti, and L. Radzihovsky, Phys. 
Rev. B 57, 7705 (1998). 

46 E.Y. Andrei, G. Deville, D.C. Glattli, F.I.B. Williams, E. 
Paris, and B. Etienne, Phys. Rev. Lett. 60, 2765 (1988). 

47 K. Moon, R. T. Scalettar, and G. Zimanyi, Phys. Rev. Lett. 
77, 2778 (1996). 

48 C.J. Olson, C. Reichhardt, and F. Nori, Phys. Rev. Lett. 
80, 2197 (1998), ibid. 81, 3757 (1998). 



22 



